Star formation in the giant HII regions of M101

The molecular components of three giant HII regions (NGC 5461, NGC 5462, NGC 5471) in the galaxy M101 are investigated with new observations from the James Clerk Maxwell Telescope, the NRAO 12-meter, and the Owens Valley millimeter array. Of the three HII regions, only NGC 5461 had previously been detected in CO emission. We calculate preliminary values for the molecular mass of the GMCs in NGC 5461 by assuming a CO-to-H_2 factor (X factor) and then compare these values with the virial masses. We conclude that the data in this paper demonstrate for the first time that the value of X may decrease in regions with intense star formation. The molecular mass for the association of clouds in NGC 5461 is approximately 3x10^7 Mo and is accompanied by 1-2 times as much atomic mass. The observed CO emission in NGC 5461 is an order of magnitude stronger than in NGC 5462, while it was not possible to detect molecular gas toward NGC 5471 with the JCMT. An even larger ratio of atomic to molecular gas in NGC 5471 was observed, which might be attributed to efficient conversion of molecular to atomic gas. The masses of the individual clouds in NGC 5461, which are gravitationally bound, cover a range of (2-8) x 10^5 Mo, comparable with the masses of Galactic giant molecular clouds (GMCs). Higher star forming efficiencies, and not massive clouds, appear to be the prerequisite for the formation of the large number of stars whose radiation is required to produce the giant HII regions in M101.Comment: 32 pages, 5 figures, accepted for publication in the Astrophysical Journa

Succinct Indexable Dictionaries with Applications to Encoding kk-ary Trees, Prefix Sums and Multisets

We consider the {\it indexable dictionary} problem, which consists of storing a set $S \subseteq \{0,...,m-1\}$ for some integer $m$, while supporting the operations of \Rank(x), which returns the number of elements in $S$ that are less than $x$ if $x \in S$, and -1 otherwise; and \Select(i) which returns the $i$-th smallest element in $S$. We give a data structure that supports both operations in O(1) time on the RAM model and requires ${\cal B}(n,m) + o(n) + O(\lg \lg m)$ bits to store a set of size $n$, where {\cal B}(n,m) = \ceil{\lg {m \choose n}} is the minimum number of bits required to store any $n$-element subset from a universe of size $m$. Previous dictionaries taking this space only supported (yes/no) membership queries in O(1) time. In the cell probe model we can remove the $O(\lg \lg m)$ additive term in the space bound, answering a question raised by Fich and Miltersen, and Pagh. We present extensions and applications of our indexable dictionary data structure, including: An information-theoretically optimal representation of a $k$-ary cardinal tree that supports standard operations in constant time, A representation of a multiset of size $n$ from $\{0,...,m-1\}$ in ${\cal B}(n,m+n) + o(n)$ bits that supports (appropriate generalizations of) \Rank and \Select operations in constant time, and A representation of a sequence of $n$ non-negative integers summing up to $m$ in ${\cal B}(n,m+n) + o(n)$ bits that supports prefix sum queries in constant time.Comment: Final version of SODA 2002 paper; supersedes Leicester Tech report 2002/1

A Search for Small-Scale Clumpiness in Dense Cores of Molecular Clouds

We have analyzed HCN(1-0) and CS(2-1) line profiles obtained with high signal-to-noise ratios toward distinct positions in three selected objects in order to search for small-scale structure in molecular cloud cores associated with regions of high-mass star formation. In some cases, ripples were detected in the line profiles, which could be due to the presence of a large number of unresolved small clumps in the telescope beam. The number of clumps for regions with linear scales of ~0.2-0.5 pc is determined using an analytical model and detailed calculations for a clumpy cloud model; this number varies in the range: ~2 10^4-3 10^5, depending on the source. The clump densities range from ~3 10^5-10^6 cm^{-3}, and the sizes and volume filling factors of the clumps are ~(1-3) 10^{-3} pc and ~0.03-0.12. The clumps are surrounded by inter-clump gas with densities not lower than ~(2-7) 10^4 cm^{-3}. The internal thermal energy of the gas in the model clumps is much higher than their gravitational energy. Their mean lifetimes can depend on the inter-clump collisional rates, and vary in the range ~10^4-10^5 yr. These structures are probably connected with density fluctuations due to turbulence in high-mass star-forming regions.Comment: 23 pages including 4 figures and 4 table

Far infrared mapping of three Galactic star forming regions : W3(OH), S 209 & S 187

Three Galactic star forming regions associated with W3(OH), S209 and S187 have been simultaneously mapped in two trans-IRAS far infrared (FIR) bands centered at ~ 140 and 200 micron using the TIFR 100 cm balloon borne FIR telescope. These maps show extended FIR emission with structures. The HIRES processed IRAS maps of these regions at 12, 25, 60 & 100 micron have also been presented for comparison. Point-like sources have been extracted from the longest waveband TIFR maps and searched for associations in the other five bands. The diffuse emission from these regions have been quantified, which turns out to be a significant fraction of the total emission. The spatial distribution of cold dust (T < 30 K) for two of these sources (W3(OH) & S209), has been determined reliably from the maps in TIFR bands. The dust temperature and optical depth maps show complex morphology. In general the dust around S209 has been found to be warmer than that in W3(OH) region.Comment: Accepted for publication in Journal of Astrophysics and Astronomy (20 pages including 8 figures & 3 tables

Star Formation in Massive Protoclusters in the Monoceros OB1 Dark Cloud

We present far-infrared, submillimetre, and millimetre observations of bright IRAS sources and outflows that are associated with massive CS clumps in the Monoceros OB1 Dark Cloud. Individual star-forming cores are identified within each clump. We show that combining submillimetre maps, obtained with SCUBA on the JCMT, with HIRES-processed and modelled IRAS data is a powerful technique that can be used to place better limits on individual source contributions to the far-infrared flux in clustered regions. Three previously categorized "Class I objects" are shown to consist of multiple sources in different evolutionary stages. In each case, the IRAS point source dominates the flux at 12 and 25 microns. In two cases, the IRAS point source is not evident at submillimetre wavelengths. The submillimetre sources contribute significantly to the 60 and 100 micron fluxes, dominating the flux in the 100 micron waveband. Using fluxes derived from our technique, we present the spectral energy distribution and physical parameters for an intermediate-mass Class 0 object in one of the regions. Our new CO J=2-1 outflow maps of the three regions studied indicate complex morphology suggestive of multiple driving sources. We discuss the possible implications of our results for published correlations between outflow momentum deposition rates and "source" luminosities, and for using these derived properties to estimate the ratio of mass ejection rates to mass accretion rates onto protostars.Comment: 12 pages, 11 gzipped gif figures, LaTex file and MNRAS style files, accepted by MNRAS, v2: reference typos and author affiliation have been correcte

HERSCHEL-HIFI spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.Comment: 6 pages, 4 figures, accepted by A&

The Computational Complexity of Generating Random Fractals

In this paper we examine a number of models that generate random fractals. The models are studied using the tools of computational complexity theory from the perspective of parallel computation. Diffusion limited aggregation and several widely used algorithms for equilibrating the Ising model are shown to be highly sequential; it is unlikely they can be simulated efficiently in parallel. This is in contrast to Mandelbrot percolation that can be simulated in constant parallel time. Our research helps shed light on the intrinsic complexity of these models relative to each other and to different growth processes that have been recently studied using complexity theory. In addition, the results may serve as a guide to simulation physics.Comment: 28 pages, LATEX, 8 Postscript figures available from [email protected]

The abundance of C18O and HDO in the envelope and hot core of the intermediate mass protostar NGC 7129 FIRS 2

NGC 7129 FIRS 2 is a young intermediate-mass (IM) protostar, which is associated with two energetic bipolar outflows and displays clear signs of the presence of a hot core. It has been extensively observed with ground based telescopes and within the WISH Guaranteed Time Herschel Key Program. We present new observations of the C18O 3-2 and the HDO 3_{12}-2_{21} lines towards NGC 7129 FIRS 2. Combining these observations with Herschel data and modeling their emissions, we constrain the C18O and HDO abundance profiles across the protostellar envelope. In particular, we derive the abundance of C18O and HDO in the hot core. The intensities of the C18O lines are well reproduced assuming that the C18O abundance decreases through the protostellar envelope from the outer edge towards the centre until the point where the gas and dust reach the CO evaporation temperature (~20-25 K) where the C18O is released back to the gas phase. Once the C18O is released to the gas phase, the modelled C18O abundance is found to be ~1.6x10^{-8}, which is a factor of 10 lower than the reference abundance. This result is supported by the non-detection of C18O 9-8, which proves that even in the hot core (T_k>100 K) the CO abundance must be 10 times lower than the reference value. Several scenarios are discussed to explain this C18O deficiency. One possible explanation is that during the pre-stellar and protostellar phase, the CO is removed from the grain mantles by reactions to form more complex molecules. Our HDO modeling shows that the emission of HDO 3_{12}-2_{21} line is maser and comes from the hot core (T_k>100 K). Assuming the physical structure derived by Crimier et al. (2010), we determine a HDO abundance of ~0.4 - 1x10^{-7} in the hot core of this IM protostar, similar to that found in the hot corinos NGC 1333 IRAS 2A and IRAS 16293-2422.Comment: 10 pages, 7 figure

The JCMT Gould Belt Survey: a quantitative comparison between SCUBA-2 data reduction methods

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artefacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software (starlink) but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only three times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered